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The continuous development of internet of things (IoT) technology enables 
many devices to be interconnected through the external environment. 
Meanwhile, 5G technology provides an enhanced quality of services with 
high data transmission rates, requiring IoT implementation in the 5G 
architecture. Free-space optical communication (FSO) is considered a 
promising technique that can provide high-speed communication links, so 
FSO is an optimal choice for wireless networks to fulfill the full potential of 
5G technology, providing speeds of 100 Gb/s or more. By implementing 5G 
features in IoT, IoT coverage and performance will be enhanced by using 
FSO models. Therefore, the paper proposed and investigated the multiple- 
input and multiple-output/free-space optical communication (MIMO/FSO) 
model using subcarrier quadrature amplitude modulation (SC-QAM) and 
relay stations over atmospheric turbulence channels by log-normal and 
gamma-gamma distribution under different turbulence conditions. The 
performance is examined based on the average channel capacity (ACC), 
which is expressed in terms of average spectral efficiency (ASE) parameters 
while changing the different parameters of the model. The mathematical 
formulas of ACC for atmospheric turbulence cases are calculated and 
discussed the influence of turbulence strength, the different number of relay 
stations, misalignment effects, and different MIMO configurations. 
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1. INTRODUCTION 


Free space optical communication is considered a promising technique that can provide high-speed 
communication links. Free-space optical communication (FSO) systems have always been proved necessary 
because of a license-free, cost-effective, high-bandwidth and high-security access technique [1]-[3]. There is 
a significant degradation in the performance of FSO communications caused by atmospheric turbulence 
along the laser beams propagation path [4]. Safari and Uysal [5] have derived the relay-assisted free-space 
optical communication. Kamiri and Kerari [6] has formulated and analyzed the bit error rate (BER) analysis 
of cooperative systems in free-space optical networks, with some recent studies on the channel model for 
amplify-and-forward (AF) systems [7]-[14]. 
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Further, the error performance of FSO systems using sub-carrier modulation schemes has been 
extensively investigated in [15]-[27]. On the subject of subcarrier quadrature amplitude modulation 
(SC-QAM) systems, the average symbol error probability of the single-input and single-output (SISO) FSO 
systems over atmospheric turbulence channels is presented in [17]. Recently, Hassan et al. [20] evaluated the 
sub-carrier intensity modulation in wireless optical communications. Bayaki et al. [21] investigated the 
performance of multiple-output (MIMO) FSO system under the gamma-gamma fading model using the 
modified Bessel function. Gripeos et al. [22] presented the average channel capacity (ACC) of MIMO FSO 
systems with different atmospheric turbulence conditions. However, an analyst is of the performance of an 
AF FSO system over atmospheric turbulence and misalignment fading is not yet available. 

This paper presents a theoretical analysis of the misalignment fading effects on the ACC 
performance of AF relay-assisted MIMO/FSO systems over atmospheric turbulence channels in section 1. 
Section 2 described the system model in detail. Next, the average channel capacity is presented in section 3. 
Section 4 discusses the results and concludes at the end. 


2. SYSTEM MODEL 
2.1. AF relay-assisted MIMO/FSO systems 

A typical serial multi-hop MIMO FSO system using SC-QAM signals with M transmitting lasers 
are pointing toward an N aperture receiver as depicted in Figure 1. The system comprises of source node, 
relaying node, and destination node. 
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Figure 1. The source node, relaying node and destination node of AF-MIMO/FSO systems 


All relay nodes are assumedly to receive and transmit signals concurrently in the same frequency 
band. The channel model of MIMO/FSO systems can be expressed by the MxN matrix, which is denoted by 
X = [Xinn mnei The electrical signal at the input of the QAM demodulator of the destination node can be 
described as (1) [21], 


M N cti 


re(t) = Relt)r p X | [Cmr 


m=1n=1 i=0 


+ Uny 


(1) 


where (X;)mn denotes the stationary random process of the turbulence channel from the m" laser to the n® 
PD, vay = D Vn; is the total accumulation noise. In this system model, we use an equal gain combining 


(EGC) scheme at the destination node to estimate the received signal from sub-channels, the instantaneous 
electrical SNR will be (2), 
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where c is the number of the relay station, y;,,,, is the random variable (RV) defined as the instantaneous 
electrical SNR component of the sub-channel from the m" laser to the n PD, it can be described by (3), 


0 


2 m=1 n=1 i= 


where X; denotes the stationary random process pdf the turbulence channel from the i AF to the (i+/)" AFR 
is the photodiode responsivity P, denotes the average transmitted optical power per symbol at each hop, 
k(0 < x < 1) is the modulation index, P; is the amplification power of the i AF module, Noy is the total 


noise variance at a destination node. 


2.2. Atmospheric turbulence channels and misalignment fading models 

In (2) and (3), X represents the channel state, which models the optical intensity fluctuations caused 
by atmospheric loss, atmospheric turbulence, and misalignment fading channel. With weak atmospheric 
turbulence, the probability density function (pdf) is described as (4) [23], 


2 x yet 1 In(X/X|Ay) +a 
N g ze x erfe( (X/X/Ao) ) 
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(4) 


where a = 0.5 o? + o? (° + c) and b = of (€? + c){(E* + c)}/2. 
With moderate to strong atmospheric turbulence conditions, the pdf of is described as (5) [27], 
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(5) 
where Ay = [erf(v)]? is the fraction of the collected power at radial distance 0, v is given by v = 


[rer / (,/2@,) with r and w, respectively denote the aperture radius and the beam waist at a distance z and 


č= a where the equivalent beam radius can be calculated by (6), 


Wzeq = w,(Vmerf(v)/2v x exp(—v’))? (6) 


where wz = Wo[1 + €(AL/mw3)*]*/? with wo is the transmitter beam waist radius at z = 0, £ = (1+ 
2w2)/p2 and po = (0.55C2k?L)~%/> is the coherence length, C? stands for the strength of the atmospheric 
turbulence, which is the altitude dependent and it is given by (7) [3], 
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1000 
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where h is the altitude in meters. v is the wind speed in meters per second, and C?(0) is the value of C? at 
the ground in m~3/?_ C2 varies from 1071”?m73/2 to 107'3m~*/ for weak up to strong turbulence cases, 


respectively. 


3. AVERAGE CHANNEL CAPACITY WITH MISALIGNMENT FADING EFFECTS 
Free space optical communication is wireless channel. It is a time variant and random channel. 
Therefore, to evaluate the channel capacity, we determine the average value, and is denoted (C). If the 
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channel's frequency response is known, the average channel capacity can also be expressed in terms of 
average spectral efficiency in bits/s/Hz, given by (8) [22]. 


(C) = Blog,(1 + yp) x f.(C)dI, (bit/s/Hz) m 


where T = {Ynm: n = 1,...,N,m = 1, ..., M} is the matrix of the channels, B is the channel’s bandwidth and 
fr (Œ) determined from the components Ynm.- 

We obtain the pdfs of AF relay-assisted MIMO/FSO systems over atmospheric turbulence channels. 
In weak atmospheric turbulence [23] and moderate to strong atmospheric turbulence [27] as (9), (10), 
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Substituting (9) into (8), develop the formula In(1 + y) = [2 (=1)* tT using mathematical formulas [9], 
2 
J exp(bz) erfc(az)dz = = exp (bz) erfc(az) — exp(—s)erf(— — az). In MIMO/FSO, the average capacity 
of log-normal can be represented by (11), 
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where c = a — In(VA2X?),t = 0.5lny + c. 


Substituting (10) into (8), the average capacity of the MIMO/FSO system for gamma-gamma can be 
given by (12), 
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Using function Meijer G, G2? [] , and use the function expansion In(1+ x) = Gy EW . The (12) is written 


as (13). 
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Performing the multiplication of the two Meijer G functions and applying the integral property. The ACC of 
AF relay-assisted MIMO/FSO systems with moderate to strong atmospheric turbulence can be given by (14). 
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4. NUMERICAL RESULTS 

The (11) and (14) can determine the average channel capacity for different atmospheric turbulence 
and misalignment fading. In (11) is used for log-normal distribution model in case of weak turbulence, while 
the other is used for gamma-gamma distribution model in the case moderate to strong turbulence. In the 
results, the value of average channel capacity is evaluated for value of turbulence strength (C7), different 
MIMO configurations, link distance L(L = 4000 m), different values of number relay stations c (with c = 0, 
c = 1, and c = 2), and amplification gain P,; = 2 dB. The constants and parameters are given in Table 1. 


Table 1. System parameters and constants 


Parameter Symbol Value 
Laser wavelength a 1550 nm 

Photo detector responsivity R 1 A/W 

Modulation index K 1 
Total noise variance No 107 A/Hz 
In-phase, Quadrature signal amplitudes M,, Mọ 8,4 

The number of relay stations c 0,1,2 
Index of refraction structure C 10" m? 

Link distance L 4000 m 


Figures 2 and 3 illustrate the ASE of average SNR with respect to y, for three values of the 
misalignment displacement standard o,, with link distances L = 4000 m, for different channel configurations 
and amplification gain P4p = 2 dB in the case moderate to strong, C2 = 3 x 10-14m~?/%. It is observed that 
the average channel capacity depends on the channel configurations and misalignment displacement 
standard, especially when increasing channel configuration and high regions of SNR. The ASE under the 
turbulence conditions is higher than in the cases of moderate and strong turbulence. It has been observed that 
with an increase in the value of c, (Figure 2, c=10, and Figure 3, c=1) capacity performance of the system 
deteriorates. 
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Figure 2. ASE versus average SNR for 8x4 QAM 
with different MIMO configurations, SISO, 2x2 and 
4x4, misalignment displacement standard o, and 
c=0, L=4000 m 


Figure 3. ASE versus average SNR for 8x4 QAM 
with different MIMO configurations, SISO, 2x2 and 
4x4, misalignment displacement standard o, and 
c=1, L=4000 m 
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Figure 4 illustrate the ASE of average SNR with respect to y, for three values of the misalignment 
displacement standard o,, with link distances L= 4000m, for different channel configurations, 
amplification gain Pp = 2dB, number relay station c=2 for strong atmospheric turbulence, 
C2 = 3 x 10714m? The results in this case change significantly from Figure 4 when the number of relay 
stations is increased by c=2. Specifically, as the value of the relay station increases, the value of the channel 
capacity decreases, this change is larger for the higher configuration. Figure 5 illustrates the ASE versus the 
number relay station in the case weak to strong atmospheric turbulence and channel configurations with 
distance L= 4000 m. The results are in agreement with the conclusions in the figures, it represents the 
dependence of the channel capacity on the system parameters. 
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Figure 4. ASE versus average SNR for 8x4 QAM Figure 5. ASE versus the number relay station, cwith 
with different MIMO configurations, SISO, 2x2 and different MIMO configurations, SISO, 2x2 and 4x4, 
4x4, misalignment displacement standard o, and link distance, L=4000 m, operational wavelength 

c=2, L=4000 m 1500 nmand 8x4 QAM scheme 


5. CONCLUSION 

This paper has theoretically analyzed misalignment fading effects on the average channel capacity 
performance of AF relay-assisted MIMO/FSO systems. The gamma-gamma and log-normal distribution 
models were used to describe the optical propagation fluctuation over atmospheric turbulence channels. The 
MIMO/FSO model is examined using subcarrier quadrature amplitude modulation and relay stations over 
atmospheric turbulence channels. The article has built an analytical expression of channel capacity for the 
cases of log-normal and gamma-gamma for MIMO/FSO system. Based on the theoretical expressions for 
ASE, the performance of SISO and MIMO systems is strongly affected by the number of AF relay stations, 
MIMO configurations, and turbulence conditions. 
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